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1 Introduction 



The low energy dynamics of strong interactions is usually investigated in the framework of 
effective theories introduced phenomenologically Now that the strong interaction is be- 
lieved to be described by quantum chromodynamics (QCD), these effective theories, if they 
are correct, should be derived from QCD in some way or another. For example, although the 
chiral perturbation theory or the Skyrme model describes the low energy physics involving the 
Nambu-Goldstone (NG) bosons, we cannot determine definitely the forms of interactions or 
the magnitude of the parameters. Then a natural question arises of how to obtain directly 
such a hadronic effective theory from QCD, which may be answered after solving the problem 
of confinement and dynamical chiral symmetry breaking. For the realistic case, QCD3 + i, this 
problem is still far beyond our knowledge. Two-dimensional QCD, however, may serve as a 
tractable model at least to know some consequences after confinement. 

Two dimensional QCD occupies a very unique position as a toy model for the 4 dimensional 
QCD. Originally it was introduced by 't Hooft 0] to demonstrate that we can sum up the 
planar diagrams and can get nonperturbative results in the limit of a large number of colors 
(the large N limit). The Bethe-Salpeter equation for the quark-antiquark system becomes 
simple and leads to what is now called the 't Hooft equation which determines the mesonic 
spectra in the large N limit. Second, two dimensional gauge theories have a linear potential and, 
thus, are confining in the lowest perturbation theory. Also, nonperturbatively, 't Hooft showed 
that the mesonic spectra are discrete, which suggests the confinement of quarks in the large 
N limit. Third, Kikkawa || presented the idea that the Hamiltonian can be easily rewritten 
only in terms of a gauge-invariant bilocal operator due to the absence of propagating degrees 
of freedom in two dimensional gauge theories. This operator is a path-ordered product with 
its ends attached by a quark and an anti-quark, and thus physically close to a mesonic state. 
Several works using gauge-invariant bilocal operators have been performed in two dimensional 
QCD ||, [5], |, [7], [|. From these facts, two dimensional QCD can be understood as an 
appropriate model to investigate how hadronic theory can be constructed after confinement 
and to discuss the higher order effects of the 1/N expansion, which is necessary for the realistic 
case N = 3. 

Since both of the works above greatly benefited from the light-front (LF) frame, we also 
work on the light-front to go beyond them. Recently the light-front Hamiltonian field the- 
ory has acquired renewed interest with the hope of solving nonperturbatively the relativistic 
bound-state problems [pll 111], |12l [13], |I4"||. Its main strategy is the "trivial" vacuum and the 



approximation to a restricted Fock space with a few number of particles. Bound states are 



2 



analyzed within a subspace with a few numbers of particles. In this sense, this approach is 
similar to the constituent quark picture, where the dynamical variables are a few constituent 
quarks. On the LF, a careful treatment of the zero mode is important for an understanding 
of spontaneous symmetry breaking (for example, see |T5|, [HJ) and vacuum structures such as 
the 9 vacuum, but we do not consider the zero modes in this paper and we will work on the 
infinite light-front space. This is reasonable because in two dimensions there is no spontaneous 



symmetry breaking for finite N theories |17[] and also there is no topological configuration in 
SU (JV) or U (N) QCD 1+1 with fundamental fermions. According to the calculations of the 
various two-dimensional models |18|, [19|, ^0], |2T, |23|, [24]] , we know that the Tamm-Dancoff 



approximation on the light-front (Light-Front Tamm-Dancoff approach, LFTD) |25], ^6], ^7|| is 
very good for lower excitations. Especially the lowest meson state is found to be a two-body 
state (a quark-antiquark state) almost perfectly. Although the success of LFTD will be partially 
due to the triviality of the vacuum, why LFTD is so good is not thoroughly understood yet and 
it is not certain that we can also describe the NG bosons in this approximation. Indeed, in the 
constituent quark model, pions cannot be described by few constituents. Since there are no NG 
bosons in two dimensions, the success of LFTD in two dimensions does not tell anything about 
the NG bosons. Furthermore, a more fundamental question of whether we can describe the NG 



bosons on the LF is not so clearly answered to date ||28|| . Therefore approaches independent of 
the Tamm-Dancoff approximation will be needed to understand why LFTD is good and to go 
to higher dimensions where the NG mode may exist. 

Motivated by these we introduce a method to obtain a bosonic Hamiltonian from fermionic 
systems. This is a well-known technique in quantum many-body physics and is called the 
boson expansion method or the boson mapping P5| . Using this method we can introduce 
bosonic variables to represent bifermion operators. On the other hand, the work of Kikkawa 
can be understood as the operator theory of the quark-antiquark bound states. The variable in 
the formalism is the bilocal operators. We extract bosonic variables from the bilocal operators 
and construct a bosonic theory. Since the main purpose of this paper is to indicate that the 
use of the boson expansion method is very natural to obtain a bosonic theory on the LF, its 
conceptual aspects are emphasized throughout the paper. 

This paper is organized as follows. In the next section, the basics of QCD 1+1 on the LF is 
presented and the gauge-invariant approach by Kikkawa is given. It is shown that an equation 
of motion for the two-body state becomes the 't Hooft equation in the large N limit. In Sec. 
Ill, we consider the general structures of the two-body states and why we need boson expansion 
methods is explained. In Sec. IV, the boson expansion methods are shortly reviewed and applied 
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to 2D QCD. Here we introduce bilocal boson operators to represent fermion bilinear operators. 
Local boson operators are given in Sec. V as collective bosons of the bilocal operators. Finally 
Sec. VI is devoted to a summary and discussion. The Appendix shows some relations for 
operators between the gauge-invariant approach and the usual one. 

2 Light-front QCDi + i and its bilocal formulation 

In this section, we fix the notation for two dimensional QCD on the light front and introduce 
its bilocal formulation of Kikkawa |Q. This formalism gives a Hamiltonian only in terms of 
gauge-invariant bilocal operators, which are essentially quark bilinears. In the large N limit, 
the Heisenberg equation of motion for the bilinear operator becomes the 't Hooft equation [0]. 

2.1 QCDi + i in light-front coordinates 

The model we consider is a two-dimensional gauge theory with a massive quark in the funda- 
mental representation of SU(N) or U(N). For simplicity we treat a one-flavor quark, but the 
generalization to many flavors is straightforward. The Lagrangian density is 

£ = -Itr(F^F^) + *(i 7 ^ - m)% (2.1) 

where F^ v = d^A u — d u A^ — ig[A^, A u ] is the field strength and = <9 M — igA^ = <9 M — igA^T a 
is a covariant derivative. The N x N matrices T a , generators of the group, are normalized as 
ti(T a T b ) = ^5 ab and satisfy the Lie algebra [T a ,T b ] = if abc T c . The fundamental quark \1/ is a 
two-component spinor and we define 

= 2 ~ k {t)' (2 - 2) 
where i — 1, . . . , N is a color index. We use the following 7- matrices; 70 = 07, 7i = ^2- 

We work in the light-front coordinates, x ± = (x° ± x 1 )/a/2, and thus regard x + as 'time' 
and so d + = gf + (= d~) as 'time' derivative. To avoid confusion, we shall use only d + and the 
spatial derivative <9_ instead of d~ or d + . If we choose the light-front 'temporal' gauge, A + = 0, 
the Lagrangian becomes 

£ = tr(<9_A + ) 2 + i^D + ij + z X f £L X - ~^{^X + XH). (2-3) 

It is evident that A + and x are n °t dynamical, because there is no kinetic term for them in 
eq. (|2.3|) . The constraints, 

dlA a + - g^T a ip = 0, (2.4) 

ill 

id_ X --y=^ = 0, (2.5) 



are solved easily and the solutions are 



2 J- 



oo 



A%{x) = 9 -j dy\x-y\J a (y)+C a x + D a , (2.6) 



oo 



1TI f°° 

where J a (x) = ip^ {x)T a %l){x) and e(x) is a sign function e(x) = x/\x\.f\ As was pointed out 
by Bars and Green ||30|| , we can set C a = D a = in the physical space. And also, since we 
do not consider the fermionic zero mode, we did not write the boundary terms in eq. (|2.7|) . 
Inserting these solutions into the Hamiltonian and using the same trick for the product of the 
sign functions as that of Kikkawa || , we have 

2 

H = J dxdx' \x-x'\J a {x)J a {x') 

+— J dxdx' e(x - x')ip\x)^(x'). (2.8) 

The remnant of the gauge field can be found as the linear potential in the first term. Finally, 
a quantization condition is imposed on the dynamical variable ip; 

{4>i (x) , ip] (y)} x+=y + = 5ij5(x - y) , 

{M x )^j(y)}x+= y + = {ipi(x),ipj(y)}x+= y + = o, (2.9) 

The Fock vacuum is 

4_|0 >= 4_|0 >= (Vi = 1,...,N , k- > 0), (2.10) 

where annihilation operators are defined as 

lH(x) = / -j= U e~ ik * + 4 e te } . (2.11) 
Jo \/2ti l j 



On the LF, as long as we do not consider the zero mode, the Fock vacuum |0 > is also the 
vacuum of the interacting Hamiltonian ( |2.8j ). 



2.2 Bilocal formulation— Gauge invariant variables 

Following Kikkawa ||, let us introduce a color-singlet bilocal field and its Fourier transform; 



1 N 

M(x,y) = -Y\^{x)My)l (2-12) 



2 1 i 



2 From now on we write the light-front spatial variables without their suffices, i.e. x and p instead of x or 
p_. Moreover, since we work in the Hamiltonian formalism, we omit the light-front 'time' coordinate. 
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/oo rjqr, rlnrrt 
-j=-j= M(x u x 2 )e^ + ^. (2.13) 
-oo y ZTT V Z7T 

Although this is defined as the commutator of fermion fields, it is not essential because we treat 
only normal ordered operators in the following discussions. If we choose the A' = gauge as 
in Ref. [[|, M is given by 

1 



M{x,y) o 



(2.14) 



where P denotes the path-ordered product. Since we are working at a fixed time in two 
dimensions, there is no path dynamics for the Wilson line. This operator is gauge equivalent 
to Eq. fl2.12|) . Indeed if we perform the gauge transformation, U{x) = P exp{ig f~°° A^dz~}, 



Eq.(^JJ) is transformed into Eq.( g7g) . 

Quantization of M(pi,p 2 ) = M*(p 2 ,pi) is given by the commutation relation; 

[: M(p l ,p 2 ):,:M(q l ,q 2 ) :] 

= :M(p 1 ,q 2 ):6(p 2 + q 1 )-:M{qi,p 2 ):5(p 1 + q 2 ) (2.15) 
+ N5( Pl + q 2 )5( qi + p 2 ){9(p 1 )9(p 2 )9(-q 1 )9(-q 2 ) - 9{-p 1 )9{-p 2 )9{q l )9{q 2 )}. 



This is derived from the anti-commutation relations ( [2.9|) . The normal order of M is defined 
through the quark operators ( |2.11 ) in terms of the Fock vacuum. The last c-number term that 



emerged from normal-ordering is very important for later discussions. Decomposing M into 
four parts by the sign of the momentum as 

M TlT2 { Pl ,p 2 ) = 9(r lP i)9(T 2 p 2 )M(p 1 ,p 2 ), (2.16) 

we can see that only : M_„ : does not annihilate the vacuum; 

: M ++ {p u p 2 ) : |0 >=: M + _( Pl ,p 2 ) : |0 >=: M_+(p 1 ,p 2 ) : |0 >= 0. (2.17) 

The state : M : |0 > corresponds to the quark- ant iquark two-body state. 3 Although the full 

Fock space is constructed by further introducing gauge-invariant baryon operators ||, we treat, 
in this paper, only the bosonic sector ( i.e. a subspace with a zero fermion number ); 



^ B =|n(:M__(p i ,g i )0|O>|. (2.18) 



The Hamiltonian fl2.8p is rewritten by using the normal ordered bilocal fields as 



Q 2 f°° 

H = -— / dp l dp 2 dq 1 dq 2 S(px + p 2 + q 1 + q 2 ) x 

07T J-oo 



3 We call this state as a "two-body state" instead of a "mesonic state" as in Ref. ||, because this is not an 
eigenstate of the Hamiltonian as we shall see. 
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P + V so ] - M (PuP2) ■■M(q 1 ,q 2 ) : 



UP1 + 92) 2 N(p 1 +p 2 ) 2 

where V denotes the principal value prescription of the integral to avoid infrared divergences 
0. The term with 1/N in the first term is absent for the U(N) theories. This Hamiltonian is 

not completely normal ordered, e.g. : M + _ :: M : is not in the normal order. This is because 

we want to keep in touch with the unit : M : which is close to the mesonic object. Although 
the Hamiltonian (|2.19|) is quadratic in terms of : M :, the Heisenberg equation of motion for 

the two-body state : M : |0 > has interaction terms due to the nontrivial algebra (|2.15|) . For 

U (N) theory (hereafter we treat only U(N) case for simplicity) we obtain 

id + : M__(g l5 g 2 ) : |0 >= [: M__(g l5 g 2 ) :, : H :]|0 > 



9 2 N r 
4tt J 


V 
dk— : 


\{ m2 


9 2 N 


2vr 



- + -] : M__(gi,g 2 ) : |0 > 
,9i 92/ 



^ J dkdp^(:M—(q 1} k-p) :: M_„(p,q 2 -k) : 



- : M__(p,g 2 ) :: M__( gi - k, k - p) : )|0 > . (2.20) 



Since it is reasonable to treat : M±± : and : M± T : to be Oi^N 1 / 2 ) and O(N ) respectively, 
we find that the first two terms give the leading contribution (0(N°)) in the large iV limit 
(N — > 00 with c/ 2 A^ fixed), and that the last term is in the next order. The last term represents 
the dissociation of a single two-body state into two two-body states. The leading contribution 
becomes a homogeneous equation and leads to the 't Hooft equation which determines the 
mesonic spectra in the large N limit, 

where /i 2 = 2r_r + is the invariant mass, r_ = q\ + g 2 is the total momentum, z = gi/r_ is the 
momentum fraction of a quark, and =< 0|(: M__(g 1; g 2 ) > is the wave function of 
the relative motion. Therefore the two-body state can be identified with a meson in the large 
iV limit. The meson does not dissociate into two mesons in leading order. 

3 General consideration of two-body states 

In the previous section, we learned from Kikkawa's formalism that even though we start from 



the Hamiltonian ( |2.19|) quadratic in terms of the bilocal operators, the resulting Heisenberg 
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equation of motion (|2.20| ) shows that there exists an interaction. In the large N limit, however, 



the interaction vanishes and thus the two-body state becomes free and is identified with the 
meson. The origin of the difference between finite and infinite N theories is that while the 
algebra ( |2.15p is nontrivial for finite N, it goes to a bosonic commutator in the large N limit 



(: M ++ : is the hermitian conjugate to : M :, see the Appendix); 

lim [: M ++ (p h p 2 ) :, : M__( ?1)?2 ) :] = N5( Pl + q 2 )5(p 2 + Ql ), (3.1) 

thus the quadratic Hamiltonian gives a homogeneous equation of motion. Therefore we insist 
that the interaction between two-body states has been hidden in the nontrivial algebra of the 
bilocal operators. What we want to do is to construct an interacting meson theory for finite 
N . It would be more convenient if we can see the interaction manifest in the Hamiltonian. 
This will be achieved if we can express the two-body operators in terms of, say, some bosonic 
operators. Using bosonic operators is justified as follows. 

For finite N theories, the real mesonic state may contain, in general, many-body components 
in addition to the two-body state; 



|meson >= \qq > +\qqqq > +\qqg > H — • , (3.2) 



where \qq > is a quark- ant iquark two-body state, \qqqq > is a four-body state, and so on. As 
was mentioned in the Introduction, in the context of the light-front Tamm-Dancoff approach 



|25| , |26| , |27| , however, we know that in various models ||18||-||24j| the lowest meson consists of 
only a two-body component. Thus to regard the two-body state as a meson would not be 
so bad at least for the lowest mesonic excitation even for finite N theory. A general (color 
singlet) two-body state with its total momentum P may be given as the linear combination of 
: M__(p, q) : |0 > with some normalization factor A/", 

\qq> = f dk-(p(k-)<$_<?}_ k _\0 > 

i J° 

= M [ P dk-<f>(k-) : M (-£;_, -P + : |0 >, (3.3) 

Jo 

where <p(k) is the momentum space wave function of the relative motion between a quark and an 
antiquark (see the Appendix for the relation between M and the quark or antiquark operators). 
Even if this state corresponds to the lowest meson, it is not a bosonic state because of the 
nontrivial algebra ( |2.15| ). We cannot construct bosonic operators only from quark-antiquark 
operators. For massless U(l) gauge theory, a regularized current operator satisfies bosonic 
commutators and it gives the usual bosonization; we obtain a Hamiltonian of a free massive 
boson. However, it contains a) a and d)d terms and is not a pure quark-antiquark operator. 
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The physical meaning of the non-bosonic properties of the two-body state is given as follows. 
Since the two-body state consists of fermions, the Pauli principle works and forbids the two- 
body states to be in the same state, for example the same momentum state. This effect works 
between two-body states. Indeed a single two-body state can be considered as a bosonic state, 
which is observed as 

< 0|[: M ++ ( Pl ,p 2 ) :, : M__( gi) q 2 ) :]|0 >= N5( Pl + q 2 )5(p 2 + (3.4) 

for arbitrary N. Hence the interaction between two-body states can be understood as the effect 
of the Pauli principle and thus this is a universal property of the fermion-pair systems. Note 
that the equation (|3.1|) can also be understood in the context of the Pauli principle; when 
iV — > oo, the number of allowed states for a fixed momentum goes to infinity and thus the 
Pauli principle loses its efficacy. 

Since the two-body state is a good approximation for the lowest meson, we want to regard 
the state as a bosonic state. Indeed it behaves as a boson in the large N limit ( [3.1|) and even 
for finite N, as long as there is a single two-body state (|3.4|) . Therefore it is natural to try to 
find a representation where : M : is expressed by some bosonic operator which is equivalent to 
the two-body state in these cases. If we are able to do this, we obtain a bosonic theory with 
interactions manifest in the Hamiltonian. As for the problem of how to represent operators 
whose algebra is nontrivial, in terms of bosonic operators, there is an answer as a standard 
technique in the many-body physics. It is the boson expansion method. Using this method, we 
can express the bilocal operators : M : in terms of bosonic operators which correspond to : M : 
in the large N limit. This method is introduced in the next section and applied to the algebra 
of : M :. 



4 Application of boson expansion method 



The boson expansion method |29| is one of the traditional techniques in non-relativistic many- 
body problems. In various systems we encounter a situation that there are some boson-like 
excitations even if the fundamental variables of the system are not bosons. The boson expansion 
method has been introduced for describing such bosonic excitations in non-bosonic systems. In 
this section we briefly review this method and apply it to our case, QCDi + i. 

4.1 Boson expansion method — a short review 



Originally the boson expansion was invented by Holstein and Primakoff pl[ to describe in- 



teracting spin waves in the Heisenberg ferromagnet. Although the first excitation from the 
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spin-aligned vacuum, i.e. a spin wave, is solved exactly, the next excitation is difficult to obtain 
due to the interaction between spin waves, which originates from the SU(2) spin algebra. So 
they introduced an boson operator and represented the spin operators by it so that they satisfy 
the 577(2) algebra. Then the theory is translated into a bosonic theory with interactions. A 
single boson state corresponds to the single spin wave. Later another representation of the spin 
algebra was proposed by Dyson ||32|| . Also, Usui |33| applied this method to plasma oscillations. 

In nuclear physics, since bosonic excitations from the ground state are very important for 
understanding the collective vibrations, the boson expansion method has been investigated in 
great detail |29], ^ |35|, |36| . It would be particle-hole operators that play an important role for 
describing (collective) bosonic excitations in the fermionic many-body problem. Let aj, be the 
creation operator of a particle with its quantum number \i and d\ that of a hole in state i. The 
four bifermion operators form an algebra djCi^, atdj, d\dj and o)^a u ; 

[didfj,, aldj] = S^Sij — Sf^djdi- Sijala^, (4.1) 

[dia^,ala a ] = 5^ v dia a) (4.2) 

[a^dv, a\a T ] = 8 ucr dla T - S^a^au, (4.3) 

[did^dldk] = SijdkCLp, (4.4) 

[d\dj,d\di\ = Sjkdjdi - Sadldi. (4.5) 

Their first motivation is to obtain collective bosons for describing the collective excitation. The 
term "collective excitation" means that the number of the participating particles is large. Thus 
they usually start with collective fermion-pair operators such as W = Ylui Cvi a td>i , and then 
represent these by bosonic operators [0. Here however we do not follow the same path for 



reasons mentioned later and here we introduce the bosonic expressions for pure fermion-pair 
operators. 

There are several ways to implement this algebra. Here we list three well-known represen- 
tations. The different boson expansions below have been shown to be equivalent [07 1- Let us 



begin with the Holstein-Primakoff (HP) type expansion. Introducing the boson operators, 

[B,*, Btj] = S^Sy, [B^ B uj ] = 0, (4.6) 
the four operators are represented by 

a t a - = J2 b Ij b ^ = a ^ ( 4 - 7 ) 

3 

Adi = Y, B lkB,u (4.8) 
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d i a » = ZXVl - A)»v B vii ( 4 -9) 

V 

= E^fv^)^ (4-10) 

V 

These equations should be considered as follows; if we substitute the r.h.s. into the algebra 
Q4.1|) - (|4.5D ) they hold by using the bosonic commutators ( |4.6|) . As long as we stay in the zero 
fermion-number sector, this replacement is correct. The square roots in (|4.9|) and (|4.10|) are 
defined by the corresponding Taylor series (*4.„„ = Y^pA^pApv)- The term "expansion" comes 
from the fact that the HP type contains an infinite expansion in terms of the boson operator. 
There exists a subspace of full bosonic Fock space, where the square roots are well defined. 
This space is called the physical subspace. If we truncate the infinite series by a few terms, it 
would not be a good approximation as it is. This is why they consider the collective state first 
and apply the boson expansion to it instead of the pure bifermion operators. 

Second, there is a finite expansion method, the Dyson expansion. In this method, the 
expressions (|4.9|) and ( |4.10| ) are replaced by o^a^ = B^ and aXd\ = B li(^~ ^)L> respectively. 
While the expansion is finite, it is evident that the relation {dia^ = aid} does not hold. 
To remedy this, the Dyson expansion needs a particular method, which is called the unitary 
projection. 

The third representation is the Schwinger type |39[, which is also finite. This method 
needs two kinds of bosons. In addition to , we introduce another boson Ay which commutes 
with B^; 

[Ay^l] = SikSji, [A i:j ,A kl \ = 0, (4.11) 
[A^B^] = [A lJ ,Bl k ]=0. (4.12) 

Then the four operators are expressed as 



ala 



= J2 B ti B ^ (4-13) 

3 

4dj = Sij-^A^A*, (4.14) 
k 

d^ = Y, A l B nk, (4.15) 

k 

a ] A = Y, B U^- (4-!6) 

k 

It is evident from Eq. ( 4.14 ) that the state annihilated by these bosonic operators is not a true 
vacuum. The true vacuum state \vac >b is given as 

\vac > B = -^Ei-^PhAhAh* ■ ■ -4^10), (4.17) 

ViV! P 
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where N is the number of the hole states, means the permutation of the indices {hi, h 2 , . . . h^}, 
and |0) is a state which satisfies B^O) = Ay|0) = 0. 



4.2 Boson expansion method — application to QCD 



1+1 



Let us go back to two-dimensional QCD. The algebra we want to simplify is Eq. (|2.15p . Essen- 
tially this is the same as that of the fermion pair operators in nuclear physics, i.e. Eqs. (|4.1| )- 
( |4.5|) . The particle-hole operators correspond to quark-antiquark operators, and so on. The 
only difference is the existence of an internal symmetry (SU(N) or U(N)) and continuum 
indices (momentum). Taking these points into account, let us introduce the bosonic operators, 

[B( Pl ,p 2 ),B^ qi ,q 2 )} = 8(p 1 -q 1 )8(p 2 -q 2 ), (4.18) 
[B( Pl ,p 2 ),B( qi ,q 2 )] = 0, (4.19) 

where momentum indices are defined to be positive (pi,qi > 0). Since : M : is constructed so 
as to be color singlet, B(pi,p 2 ) need not carry color. In this paper we consider only HP type 
expansion. The bifermion operators are expressed as 

poo 

: M. + ( Pl ,p 2 ) : = / dq B\-p l ,q)B{p 2 ,q)=A{p 2 ,-p l ), (4.20) 
Jo 

poo 

:M + 4 Pl ,p 2 ): = - dq B\q, -p 2 )B(q, Pl ), (4.21) 

Jo 

poo 

:M ++ ( Pl ,p 2 ): = / dq (VN=A)(p 2 ,q)B(q, Pl ), (4.22) 
Jo 

POO 

: M__(p 1 ,p 2 ) : = / dq B^{q,-p 2 ){VN^A){q,-pi). (4.23) 
Jo 

Since : M ++ := 0(N 1 / 2 ) and : M_ + := O(N ), the bosonic operator is considered to be 0(N°). 
Of course, the square roots in Eqs. Q4.22j ) and ( |4.23| ) are defined by their Taylor expansions, 
for example, 

:M ++ ( Pl ,p 2 ): = VNB(p 2 , Pl ) (4.24) 

dqdkB\q,k)B(p 2 ,k)B(q, Pl ) . 



2VNJ0 

As was noted in the previous section, the bosonic representation should satisfy two conditions 
to keep in touch with the meson picture. The first requirement that : M : itself becomes bosonic 
in the large N limit, i.e. Eq. Q3.lp be satisfied as 

lim -= : M__{ Pl ,p 2 ) := B\- Pl , -p 2 ). (4.25) 

N^oo </ J\j 



And also the second one that the single two-body state should be bosonic for finite N is observed 
apparently as 

: M__(p x ,p 2 ) : |0 >= VNB\- Pl , -p 2 )\0), (4.26) 
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where |0) is a vacuum state of the bosonic operator; B(p, g)|0) = 0. That is why we used the 
HP type expansion. 

Since the boson expansion becomes a good approximation for systems with large partic- 
ipating particles, they start with collective bosons instead of pure particle-hole operators in 
nuclear physics In our case, however, there is a parameter N associated with the internal 



symmetry and the bilocal operators are color singlet. The number of participating particles in 
the two-body state goes to infinity if N — > oo. Thus we do not have to introduce collective op- 
erators at this stage. Collective bosons are introduced to give local bosons in the next section. 
Furthermore, note that also in the large N limit, the effect of the physical subspace vanishes. 
The existence of internal symmetry makes the situations simple. 

Substituting ( |4.2q )- fl4T23| ) into the Hamiltonian ( |2.19| ), we obtain for the leading contribu- 



tion; 

r(o) _ 9 2 N 



(o) 9 1 » f°° 

H ( H ' P = — — / dp 1 dp 2 dq 1 dq 2 5(p 1 +p 2 - qi - q 2 ) x 

S7T JO 



v 

; B\p 1 ,p 2 )B(q 2 ,qi) 



(Pi ~ q 2 ) 2 ' 

+\ [rn 2 ~ jT dqdq' + j\ B\q, q')B(q, q'), (4.27) 

where the suffices (0) and HP imply the 0(N°) contribution and the Holstein-Primakoff type 
expansion. The equation of motion is 

id + B(p, r-p) = [B(p, r-p), Hfijp] 

g 2 N r V 



L dq T^aY B ^ r - q) 



An Jo (p — q) 

1 / a 2 N\ (V V \ 

+ - [m 2 - y — \ [- + )B(p,r-p). (4.28) 

2\ 2n J \p r — pi 



Essentially this is equivalent to the 't Hooft equation ( |2.21 



The higher order Hamiltonians {H^ P = 0(N~ n / 2 )) are obtained by substituting the expan- 
sion of Eqs. ( [4.20p -( fl.23j) . The next leading Hamiltonian is on the order of iV -1 / 2 and is given, 
after normal ordering, by 

r(D _ 9 2 Vn r^ ^ ^u ^^.^^ v 



Hrjp = — / dpxdp 2 dq\dq 2 dk 5(p\ —p 2 + qi + q 2 ) 

An Jo 



X 



Oi + Q2) 2 

{B\k,p 2 )B(k, Pl )B(q 2 , Ql ) - B\p 2 , k)B{ Pl , k)B{q x , q 2 ) + h.c.} , (4.29) 

which corresponds to a three-point vertex. 

The expressions fl4.20| )-( [4\23| ) are very similar to those of Nakamura and Odaka M. They 
introduced bosonic operators as the leading contribution of : M : in the 1/N expansion. This 
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construction is the same as the requirement Q4.25 ) and indeed their bosonic operator is equiv- 



alent to ours in the large N limit. However, if we express their representation in our notation, 
it is different from ours; 

: M_ + (p 1 ,p 2 ) : = dq B^(-p u q)B{p 2 ,q), 
Jo 

POO 

:M+_(pi,p 2 ): = - / dq B\q,-p 2 )B(q,p 1 )=C(p 2 ,-pi), 

Jo 

POO 

:M ++ (p 1 ,p 2 ): = / dq (VN + C)(p 2 ,q)B(q, Pl ), 
Jo 

POO 

:M_.(p u p 2 ): = / dq B*(q,-p 2 )(y/N + C)(q,-pi). 
Jo 

This is one of the HP type expansions. They were not aware that this is the boson expansion 
methods. However, once we recognize the usefulness of the boson expansion method, we can 
apply this idea to many other models even with small N. Other finite expansion methods such 
as the Dyson type or methods using collective fermion pairs will be found to be more useful for 
such theories. 

5 Local boson operators as collective states 

It was bilocal operators that we introduced in the previous section. Here we construct local 
boson operators. The following is essentially the same as that of Nakamura and Odaka [|]]. 
However, according to the many-body physics, their construction is understood as the intro- 
duction of the collective states from B(p, q). Collective boson operators are defined as the linear 
combination of bilocal bosons with their total momentum r, 

bt = y/f [ dz <f)(z)B\rz, (1 - z)r). (5.1) 
Jo 

Let us assume that the wave function <f)(z) forms an orthonormal complete set; 

dz 4>* n (z)(j) m (z) = 5 nm , (5.2) 
Y,UzW n {z')=5{z-z'). (5.3) 

n 

As far as {<p n } is a complete set, the operators 

&(")t = yft f 1 dz <p n (z)BUrz,(l- z)r), (5.4) 
Jo 

&(») = ^ \ dz c/)* n (z)B(rz,(l- z)r), (5.5) 
Jo 

satisfy the bosonic commutators 

[&W 6 <m)t] = SnmS{r _ r / }) [ 6 W 6 M] = o. (5.6) 
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Let us determine {4> n } so that b and W diagonalize the lowest Hamiltonian H^ p . Inserting 
these operators, the equation of motion can be written as 

( ^ 9 2 N\(l 1 \ g 2 N [i , M*') ( , 7 , 

MK{z) = [ m - — [- + — j M*) ~^rj dz T^Vv ( 5 - 7 ) 



which means that the wave function <p n {x) is a solution of the 't Hooft equation. If we define 
a scalar field by 

= -L r (feWe-*™ + fc^V™) , (5.8) 



/27T 7o V 2 ^ 

this is a field operator of the n-th excitation state. 3> n (x) satisfies the usual LF commutator; 

[$ n (a0, ^(3/)]x+=y+ = -^^ nm e(x- - y"). (5.9) 
Using this operator, H^ p is diagonalized and can be written as 

= Ef^ 4 : ®n{x? ■ ■ (5.10) 

ra J — oo Z 

This is the Hamiltonian of free scalar bosons with their mass being the eigenvalues of the 't 
Hooft equation. The next order Hamiltonian is 

g 2 VN dR /-i 

az 

nml 



ri ri Jz(l - z) 

nmi[z) = dx dy T — — — (f>* n (zx)(f> m (x)(f>i(y), 

Jo Jo {z(l — x) + (1 — z)yY 



where 

K nrr d(z) 

ri ri Jz(l-z) 

K nm i{z) = dx dy- — r- 0*(1 - z(l - x))(f) m (x)(f>i(y). 

Jo Jo \zx + (1 — z)(l — z)Y 

Hence for finite N, we obtain a Hamiltonian of infinite kinds of bosons interacting with each 
other. In principle, there are infinitely many-point vertices which originate from the infinite 
expansion of the HP type. 

Recently Barbon and Demeterfi [[| derived an effective Hamiltonian in the approximation 
that the following operators are bosonic; 

a P ee JV [ P dk-4>(k-) : M (—&;_, -P + fc_) : , (5.12) 

Jo 

a P ee M [ dk_(j>*{kJ) : M ++ (P : . (5.13) 
Jo 
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Indeed they satisfy [a P , Oq] = S(P — Q) + 0(1/ N) and thus are bosonic in the large N limit, 
which has been insisted on many times. It is characteristic that their Hamiltonian has only 
finite terms in contrast to ours. This is clearly because they did not incorporate the next 
leading order of the commutator. Their method cannot be extended to small N theories, while 
the boson expansion methods have other ways to avoid it, as has been noted. 

6 Conclusion and Discussions 

In this paper we have found that the boson expansion methods are useful for constructing 
bosonic theories from two dimensional QCD on the light-front. In this method, fermion bilinear 
operators can be nonlinearly expressed by bilocal bosonic operators such that the same algebra 
holds. Among various representations, if we choose the Holstein-Primakoff type expansion, 
the bosons are identified with mesons in the large N limit, and for finite N we obtain an 
interacting boson theory. The physical meaning of this procedure is as follows. Because of the 
inevitable effects of the fermion statistics, the two-body states cannot stay bosonic. This effect 
is called the kinematical effect because it is independent of the Hamiltonian. The kinematical 
complexity in the two-body states is translated into interactions among bosons via the boson 
expansion methods. Furthermore, we can build local boson operators from bilocal bosons as 
the collective state. Since this is determined to make the lowest Hamiltonian diagonal, this 
effect is called the dynamical effect. Eventually we have obtained a local interacting theory 
with infinite kinds of bosons. Thanks to the existence of the parameter N, we can argue the 
above two effects separately, which is not the case in nuclear physics. 

In fact, the introduction of the gauge invariant operator was not crucial for implementing 
the boson expansion method, but it has made the argument very transparent and leads to 
the discovery of the applicability of the boson expansion methods. Although there are no 
dynamical degrees of freedom in the gauge field when the space is infinite, if we work in finite 
(light-front) space, there exist dynamical zero modes of the gauge field. The gauge-invariant 
operators contain them which is not the case in nuclear physics. The boson expansion method 
can also be applied to these operators as far as the winding of the zero mode is trivial. 

The application of the boson expansion method has been performed for light-front field the- 
ory. However, we can do the same things for perturbation theory in an equal-time formulation. 
The merit on the light-front is, thanks to the triviarity of the vacuum, that the few-body states 
are expected to be good variables for describing hadrons. This is true even for small N theories. 
Thus we expect that bosons also describe hadrons. In equal-time formulations, however, we 
cannot expect in general that only the two-body state gives a meson. Thus boson expansion 
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methods are expected to work better on the light-front than on the equal-time. 

There are many things to be clarified in this formalism. Especially, it is important to 
know the relation to the usual bosonization. In bosonization, the number of bosons is finite, 
while it is infinite in our formalism. As a more concrete problem, the existence of meson- 
meson bound states is suggested in Ref. It would be interesting to investigate it in 
our framework. Furthermore, there is the problem of how we can describe baryons. Recently 
Rajeev || argued that the bilocal operator M(p, q) is an element of an infinite dimensional 
Grassmannian, and that the baryon number is understood as the topological invariant (virtual 
rank) of the Grassmannian, and thus baryons are interpreted as solitons. One of our motivations 
was to see how we can obtain an effective theory such as the Skyrme model. In the Skyrme 



model, a baryonic state is considered as a soliton in the large N limit flOfl . Our resulting 
Hamiltonian is a nonlinear theory of mesons. Although the topological effect is difficult to 
observe in the trivial vacuum of the light-front, it is intriguing to find a solitonic configuration 
in our formalism. 

In nuclear physics, the boson expansion method is considered as a method that goes beyond 
the Tamm-Dancoff approximation. Since the structure of the LF theories is similar to non- 
relativistic quantum theory, it seems to be natural that we have found that the boson expansion 
method is also applicable to LF field theory when we want to go beyond the TD approximation. 
Also, in many-body physics, the boson introduced is thought to be appropriate for describing the 
Nambu-Goldstone boson when there is symmetry breaking. Indeed at the very first application, 
the spin wave corresponds to the NG boson under spontaneous magnetization. Therefore there 
might be a possibility that also on the LF the boson introduced in the boson expansion methods 
can describe the NG boson. Since this is a very important problem, we should apply the boson 
expansion method to higher dimensional theories where the NG boson or the pseudo NG boson 
will exist. 
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Appendix 

If we use the fermion operator after solving the constraints, we can easily relate the bilocal 
operators : M : to the quark or anti-quark operators. As in eq. ( |2.11| ), the usual definition of 
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the quark or anti-quark creation/annihilation operators is given by 

poo flic f 

V*0r) = / -j= {4 e~ ikx + 4 e lkx \ , (A.l) 

■JO V27T L J 

where the momentum takes only positive values > 0. On the other hand, Fourier transfor- 
mation of the gauge invariant operators is given as in eq. ( |2.13| ). Thus its momenta can take 
negative values. Taking care of this point, the relation between two is given as follows: 

:M ++ (p,q): = d v a\ (A.2) 

:M__(p,g): = afl p <fl g (A.3) 

:M + 4p,q): = -Sl q tf p (A.4) 

:M_+(p,g): = a%a\ (A.5) 

From these we can see that : M ++ : is the hermitian conjugate to : M :, 

(:M ++ (p,g):) t =:Af__(-g,-p): . (A.6) 
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